Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 19 December 2008 (MN L3TbX style file v2.2) 



Unveiling the birth and evolution of the H ii region Sh2-173 

S. Cichowolski/* G.A. Romero,^'^ M.E. Ortega/f C.E. Cappa, 2* 
■ and J. Vasquez,^ | 

^ Instituto de Astronomi'a y Fi'sica del Espacio (lAFE), CC 67, Sue. 28, 1428 Buenos Aires, Argentina 
, ^ Instituto Argentina de Radioastronomi'a ( CCT-La Plata, CONICET), CC 5, 1894 Villa Elisa, Argentina and Facultad de Ciencias 
1 Astrondmicas y Geoffsicas, Universidad Nacional de La Plata, Paseo del Bosque s/n, 1900 La Plata, Argentina 
• ^ Departamento de Fi'sica y Astronomi'a, Facultad de Ciencias, Universidad de Valparaiso, Valparaiso, Chile 



ABSTRACT 

Based on a multiwavelength study, the interstellar medium around the H II region Sh2- 
173 has been analyzed. The ionized region is clearly detected in the optical and in the radio 
continuum images. The analysis of the H 1 data shows a region of low emissivity that has 
an excellent morphological correlation with the radio continuum emission. The H II region is 
partially bordered by a photodissociation region, which, in turn, is encircled by a molecular 
structure. The Hi and CO structures related to Sh2-173 are observed in the velocity ranges 
from -25 to -31kms^^, and from -27 to -39 kms^^, respectively. Taking into account the 
presence of noncircular motions in the Perseus spiral arm, together with previous distance es- 
timates for the region, we adopt a distance of 2.5 ± 0.5 kpc for Sh2-173. Seven hot stars were 
identified in the field of Sh2-173, being only one an O-type star. The amount of energetic 
photons emitted by this star is enough to keep the region ionized and heat the dust. Given 
that an expanding H II region may trigger star formation, a search for YSO candidates was 
made using different infrared point source catalogues. A population of 46 YSO candidates 
was identified projected onto the molecular clouds. 

On the other hand, Sh2-173 is located in a dense edge of a large (~ 5°) Hi shell, 
GSH 1 17.8+1.5-35. The possibility for Sh2-173 of being part of a hierarchical system of three 
generations is suggested. In this scenario, the large H I shell, which was probably originated 
due to the action of Cas OB5, would have triggered the formation of Sh2-173, which, in turn, 
is triggering new stars in its surrounding molecular cloud. To test this hypothesis, the ages 
of both, the H II region and the large shell, were estimated and compared. We concluded that 
Sh2-173 is a young H II region of about 0.6 - 1.0 Myr old. As for the large shell we obtained 
a dynamical age of 5 ± 1 Myr These age estimates, together with the relative location of the 
different structures, support the hypothesis that Sh2-173 is part of a hierarchical system. 

Keywords: ISM: individual:Sh2-173 - ISM: kinematics and dynamics - Hii regions - 
stars: formation. 



1 INTRODUCTION 

Massive stars affect their environment by ionizing radiation, stellar 
winds, and, finally, supernova explosions. These stars produce in- 
tense extreme ultraviolet (hu > 13.6 eV) radiation, creating Hll 
regions around them, and intense far ultraviolet (4.5 < Iw < 13.6 
eV) radiation that can penetrate into regions of molecular gas, dis- 
sociating it. They also have powerful winds which sweep up the 
surrounding gas creating what is known as an interstellar bubble, 
i.e. a minimum in the gas distribution characterized by a low vol- 
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ume dens ity and a high temp erature, surrounded by an expanding 
envelope jWeaver et aljl97'7b . In the case of OB associations, their 
collective effect has a huge impact in the surrounding interstellar 
medium (ISM), creating large structures of about hundred parsecs. 
Winds and supernova explosions transfer to the ISM vast amounts 
of mechanical energy, generating the compression of nearby molec- 
ular clouds. The effects of the massive stars are mostly destruc- 
tive in their immediate neighborhood, since they tend to disrupt the 
parental molecular cloud. However, a little further away, massive 
stars stimulate star formation. A number of papers have been de- 
voted to investigate star formation in the dense molecular clouds 
in the surroundings of H II regions, both from a theoretical and an 
observational point of view. 

Different processes can trigger stellar formation at the pe- 
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2007). Star formation 


can be induced in the environs of mas- 



sive stars due to the action of expanding shocks which mod- 
ify the surrounding material, either favouring the formation of 
cometary globules in preexistent clumps ("Radiative Driven Implo- 
sion" model, RDl) or sweeping up the surrounding material into a 
collected layer, which latter becomes gravitationally unstable gen- 
erating massive fragments ("collect and collapse" model). Both 
processes, which may act simultaneously in an Hll region, lead 
to the formation of a new generation of stars. Observational evi- 
dence of both (the RDI and "collect and c ollapse") processes has 
been found in a number of H II regions (e.g . Thompson et al .|2004 
Deharveng & Zavagnd l2008l : ICappa et aU 



20081) . 



20081 : iKirsanova et al.l 



The goal of this paper is to achieve a better understanding of 
the phenomena associated with the interaction between young mas- 
sive stars and their surrounding interstellar medium and, specifi- 
cally, to attempt to shed some light on the process of triggered star 
formation. To carry out this study a multiwavelength approach is 
essential in order to obtain a complete picture of the different com- 
ponents of the ISM, the young stellar population, and the interac- 
tions among them. 

In this work we investigate the environs of the Hll region 
Sh2-173 by analyzing the distribution of the ionized and neutral 
material, and that o f the interstellar dust. Sh2-173 is an optically 
visible Hll region jSharplesj[l959l) located at {I, b) = (119?4, - 
0?94) in the Perseus spiral arm. The DSS2-R image of Sh2-173 is 
displayed in Fig. [T] The region is approximately circular in shape 
with a diameter of about 30'. An arc-like structure of enhanced 
optical emission is evident on the north-west side of the region, 
while fainter emission is present at higher galactic longitudes and 
lower galactic latitudes. Inside the ionized region, two voids lack- 
ing optical emission are evident at (/, b) = (1 19?4, --0?85 ) and (/, b) 
= (119?55, -0?95). Russeil, A dami. & GeorgelinI l l2007h identified 
seven hot stars in the field of Sh2-173, which are listed in Tabl^ 
and indicated with small triangles in Fig.[T] All of the stars but one 
are projected close to the center of the region. ALS6155, on the 
contrary, is seen projected on its southeast border. 

The distance to Sh2-173 has been estimated by several au- 
thors using independent met hods. A spectroph otomet ric distance of 
2.7± 0.9 kpc was derived bv lGeorgelinI jl975l) , while lRusseil et all 
(12007) estimated a distance of 3.12 ± 0.34 kpc based on UBV 
photometric and spectroscopic data of the exciting stars (see Ta- 
ble [TJ. Using a new method for determination of distances, which 
is based on Hi column d ensities and is independent of Galactic 
kinematics, [Foster & Rou tledge ( 2 003h estimated for th i s H II re- 
gion a distance of 2.2 ± 0.4 kpc. lilitz. Fich, & StarkI h982l) in 
their CO survey towards H II regions found molecular gas related 
to Sh2-I73 at -34.5 kms^^ (all the velocities are referred to the 
LSR). This velocity coincides with the radial velocity of the ion- 
ized gas derived from Fabry-Perot H a observations, —34.3 ± 0.3 
kms"^ l lFich, Dahl. & Treffersll99"a) . It is well known that noncir- 
cular motions on a la rge scale are present in the Per seus spiral arm 
Brand &Bli"t3l 19931) . An inspection of Fig. 2b of iBrand & Blitzl 



19931), which takes into account noncircular motions, suggests for 
a velocity of about -34 kms^^ , a kinematical distance of 1.8-2 
kpc. Thus, the distance of Sh2-173 is between 1.8 and 3.1 kpc. As 
a working hypothesis, we will adopt a distance of 2.5 ± 0.5 kpc for 
Sh2-173. 

Sh2-173 is seen located in the dense border of a large H I shell 
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Figure 1. DSS2-R optical image of the H II region Sh2-173. The small tri- 
angles indicate the position of the seven OB stars listed in Table[T] 



Table 1. Stellar parameters of the OB stars associated with Sh2-173. 



Star 


Spectral 
Typel 


Galactic 
coordinates 


Distance^ 
[kpc] 


Comments 


ALS6145 


B0.5V 


11 9? 39,-0? 93 


3.17 ± 1.27 




ALS6150 


B21V 


119?44,-0?91 


3.78 ±1.89 




ALS6151 


09V 


11 9? 44,-0? 92 


2.57 ±0.5 


BD+60 39 


ALS6155 


B0.5V 


119?44,-1?11 


3.36 ±1.34 




ALS6156 


B2V 


11 9? 48,-0? 93 


4.41 ± 1.77 




ALS6157 


B2IV 


119?48,-1?02 


2.27 ± 1.14 


BD+60 42 


ALS6158 


BIV 


119?52,-0?85 


4.02 ±0.66 





1- From lRusseil et al.l i2007h . 



(~ 5° diam eter) that has been found i n the Perseus s piral arm by 
lFici] | |l986h using IWeaver & Williams! ( Il973l) H I data. lFichI ( 11980) 
suggested that the presence of SNR remnants and several H II re- 
gions in the fie ld indicates a great deed of s tar-formation activity go- 
ing on nearby. iKiss. Moor. & Totl? ('2004") foun d a far infrared rinj 
close to the H I shell defined by Fich (19861) . iHever et al.l ( Il99l 
found a large void of CO emission centered at (/, b) = (1 17? 6, 0?0), 
coincident with the Hi shell. Based on a multiwavelength study, 
|Mo6r & Kiss ( 20()3) determined the main parameters of the shell 
and associated i t with the OB as s ociatio n CasOB5 (= Rupr33). 
According to .Garmanv & Stencell jl992h . Cas 0B5 has a zero-age 
main sequence (ZAMS) fairly good defined with a distance modu- 
lus DM =11.5 mag (D ~ 2 kpc). The diameter of the association 
is about 4°, or 140 pc at 2 kpc. Its earliest star is HD 108 (06f). 
Among the ir members there are six 0-type stars, being the oth- 
ers B stars dGarmanv & Stencellll992l) . A WN4 star (WR 159) was 
found in the area of Cas0B5, located at (I, b) = (I15?78, +V.1A) 
( lNeguerueIall2003l) . 



Table 2. Observational parameters. 
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Synthesized besm. (0^ x 


1 ' 1 X 1 '0 


Position angle 


07° 


Velocity coverage 


-164 to 44.7 kms~l 


Channel separation 


0.82 kms-l 


Velocity resolution 


1.3 kms~^ 


Observed rms noise (single channel) 


1.3 K 


IxdUlU L.UII11I1UUII1 dl H'Uo 1V1X].£. 




Synthesized beam (0^ x ©jj) 


3'2 X 2'8 


Position angle 


95° 


Observed rms noise 


0.5 K 


KaQlO COniinUUln al ivlxlZ. 




Synthesized beam (0^ x ©j, ) 


yl. yj X u. o 


Position angle'* 


y / 




0.05 K 


CO data: 




Angular resolution 


1.'67 


Velocity coverage 


-150.7 to 38.9 km s-i 


Velocity resolution 


0.98 kms-i 


Channel separation 


0.82 kms-i 


Observed rms noise (single channel) 


0.2 K 


Infrared data: 




HIRES 




Angular resolution 


0f5 - 2' 


MSX 




Angular resolution 


18^4 



a Position angle measured CCW from the longitude axis. 

2 DATA SETS 

Radio continuum data at 408 and 1420 MHz, as well as 21 -cm 
spectral line data, were obtained using the Dominion Radio Astro- 
physical Observatory (DRAO) interferometer as part of the Cana- 
dian Galactic Plane Survey (CGPS) ( Taylor et al. 2003) . The COPS 
data base also comprises other data sets that have been reprojected 
and regridded to match the DRAO images. Amon g them ar e the 
IRAS high resolution (HIRES) data (Fowler & Au mannlll994h and 
the Five College Radio Astronomi cal Observ atory (FCRAO) CO 
Survey of the Outer Galaxy (Hever et al.lll998h . 

Radio con t inuum surveys at 325 MHz (W ENSS) 
llRengelink et al.' 'l997^, 2695 MHz dFurst et al.l Il990l) . and 
4850 MHz (Condon et al. 1994) were also used in this work. 

Mid -infrared data w ere taken from the MSX Galactic Plane 
Survey JPrice et al. I I2OOII) and near-infrared data were obtained 
from the 2MASS All-Sky Point Source Catalog jSkrutskie etaTI 
l2006h . 

Table |2] summarizes the most relevant observational parame- 
ters. 

3 SH2-173 AND ITS LOCAL ISM 
3.1 Ionized Gas 

Figure|2]shows the radio continuum emission distribution at differ- 
ent wavelengths. The upper panel shows the radio emission at 408 
MHz, while the image at 1420 MHz is shown in the middle panel. 
Clearly, the radio emission at 408 MHz is dominated by strong 
point like sources. An arc of strong radio continuum emission over- 
laying a more extended diffuse emission is observed at 1420 MHz. 
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Table 3. Radio continuum flux densities of Sh2-173. 



Frequency (MHz) 


Integrated flux density (Jy) 


325 


1.8 ± 1.0 


408 


1.3 ±0.6 


1420 


1.8 ±0.3 


2695 


1.8 ±0.3 


4850 


1.6 ±0.3 



The radio emission extends towards higher galactic longitudes and 
lower galactic latitudes, slowly decreasing in brightness. The bot- 
tom panel displays an overlay of the emission at 1420 MHz (line 
contours) and the optical image (grey-scale). The arc-like structure 
is coincident with the brightest optical region. On the contrary, the 
relatively strong radio emission near (I, h) = (119°4, -0?85) is pro- 
jected onto one of the regions lacking optical emission. The H II 
region shows a similar morphology at 325, 2695, and 4850 MHz. 

To elucidate the nature of the numerous point sources ob- 
served in the field, we have used the av ailable catalogued flu x 
densities (The Texas Survev [365 MHz. JPouglas et alj|l996h l. 
WENSS [325 MHz jRengeUnk et"ai]| 19971) 1. NVSS [1420 MHz, 
( ICondonetal.l ll998')1). The derived spectral indexes (Si, ~ i^") 
are representative of extragalactic nonthermal radio sources. In par- 
ticular, the strong radio source located at (/, h) = (119? 38, -0?84), 
which is close to the bright part of Sh2-173, has a flux density of 
235 mjy at 365 MHz and of 67.1 mJy at 1420 MHz, resulting in a 
spectral index a ~ —0.9. The contribution of the compact sources 
projected onto the H II region was subtracted in order to properly 
estimate the flux density of the ionized region. The derived flux 
densities of Sh2-173 are listed in Table [3] The least square fitting 
of these values gives a spectral index a = 0.0 ± 0.1, confirming 
that the extended source is an H II region. 

The parameters of the ionized gas were derived from the im- 
age at 1420 MHz. The electron density and the H II mass we re ob- 
tained using the expressions of Mezg er & HendersonI l ll967^ for a 
spherical H II region of constant electron density (rm s electron den- 
sity H e), adopting an electron temperature of 9000 K dOuireza et al.l 
l2006h . The derived parameters are listed in Table|4] 

The number of UV photons necessary to keep the gas 
ionized was derived using NLy-dW'''' s-^) = 3.51 X 
10~^ 7ie(cm~^) Rg (pc). Errors in the linear radius, in the rms 
electron density, and in the number of UV photons come mainly 
from the distance uncertainty. 

Bearing in mind that massive stars have a copious UV 
flux capable of ionizing the surrounding Hi gas, we investi- 
gate whether the massive stars in Sh2-173 can provide the en- 
ergy to ionize the gas. As described in Sect. [T] the only O- 
type star among the exciting stars is BD-l-60 39, being the other 
stars B ones. Given that the number of ionizing photons emit- 
ted by an 0-type star is significantly larger than that of a B- 
type star, we neglect the contribution of the UV photons emit- 
ted by the stars ALS 6145, 6150, 6155, 6156, 6157, and 6158. 
The UV photon flux correspondi ng to an 09V star is in the range 
N-^ (s~^) = (1.1 - 3.6) xlO^*^ (IVacca. Garmanv. & Shun 199^; 
iMartins. Schaerer. & Hillieil200^ . By comparing A^* with the UV 
photons needed to keep the gas ionized (NLy-c) listed in Table|4] 
we conclude that the 0-type star alone can maintain the H II region 
ionized. 
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Table 4. Main parameters of the observed structures related to Sh2-173. 
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Figure 2. Radio continuum images of the Hll region Sh2-173 obtained 
from the CGPS. The small triangles indicate the position of the OB stars, 
(a) 408 MHz. Contours levels are at 58 and 62 K. (b) 1420 MHz. Contour 
levels are at 5.7 and 6.8 K. (c) DSS2-R image (grey-scale) with the 1420 
MHz radio continuum contours superimposed. The plotted area is smaller 
than the one shown in the other images. 



Adopted distance (kpc) 


2.5 ± 0.5 


Ionized gas 




Flux density at 1420 MHz (Jy) 


1.8 ±0.3 


Angular radius (arcmin) 


13 


Linear radius (pc) 


9.5 ± 2.0 


Ionized mass ) 


480 + 250 


rms electron density rie (cm~^) 


6± 1 


Emission measure (pc cm~^) 


1000 ± 50 


Used Lyman UV photons [N^y-c] {^~^) 


(9.3 ±4.0) xlO""" 


IR emission and dust parameters 




Si2 (Jy) 


50 ± 12 


S25 (Jy) 


67 ± 30 


Seo (Jy) 


590 ± 170 


Sioo (Jy) 


1 100 ± 500 


^IR (L©) 


(9±4)xl03 




J 7 -|_ Q 7 


Td (K) 


30 ±5 


H I structure 




Central coordinates 


119? 45, -0?93 


Angular radius (arcmin) 


14 


Linear radius (pc) 


10 ±2 


Velocity interval (vi,V2, kms~^) 


-31.1, -24.6 


Systemic velocity (kms^^) 


-27.0 ±1.3 


Expansion velocity (kms^^) 


5 ± 2 


Velocity width (Ad, kms^^) 


6±2 


Neutral atomic mass in the shell, Mgh (Mq) 


150 ± 100 


Missing atomic mass in the cavity, Mmiss (Mq) 


120 ± 80 


Molecular clouds 

Mean temperature (K) 
Mean area (pc^ ) 

Velocity interval (vi,V2, km s~ ^ ) — 

Systemic velocity (kms~^) 

Mean H2 column density (10^^ molcm"^) 

Molecular mass (10^ Mq) 

Volume density riu (cm""^) 


Cloud A Cloud B 
0.9±0.1 1.2 ±0.1 

80 20 
38.6, -27.8 -39.4, -34.4 
-32.4 ±1.0 -36.5 ±1.0 
1.9 ±0.6 1.1 ±0.3 
5.8 ±2.0 1.0 ±0.3 
520 ± 210 560 ± 175 



3.2 The Dust 

The HIRES IRAS data have been searched for heated dust counter- 
part to the ionized hydrogen region. Figure[3]shows the 12, 25, 60, 
and 100 /im images. The image at 12 /im reveals a circular struc- 
ture which exhibits less emissivity in its central part and is brightest 
towards lower galactic longitudes. A more diffuse emission distri- 
bution is observed at 25 /xm. The brightest side at this wavelength 
coincides with the strongest radio continuum emission. The images 
at 60 and 100 /^m exhibit morphological features similar to those 
observed at 25 /^m. The observed IR morphology and the position 
of the OB stars suggest that the dust is being heated by the stellar 
radiation. 

The left panel of Fig.|4]exhibits the mid-IR distribution at 8.3 
^m (MSX Band A). The image resembles the IR distribution ob- 
served at 12 /xm. In the right panel a comparison between the IR 
emission at 8.3 /jm (grey-scale) and the radio continuum emis- 
sion at 1420 MHz (contours) is shown. Clearly, the IR emission 
is mostly seen projected outside the continuum emission. The arc 
feature observed at 8.3 ^m is not detected in the other three MSX 
bands (named C, D, and E). These facts suggest that the polycyclic 
aromatic hydrocarbons (PAHs) could be the main responsible for 
the emission detected at 8.3 /im, indicating the existence of a pho- 



todissociated region (PDR) in the border of the Hll region. The 
absence of 8.3 /xm emission in the interior of the ionized region 
can be interpreted as the destruction of PAH molecules by the ex- 
treme ultraviolet radiation of the ionizing stars. The fact that there 
is emission detected in the 25 fim IRAS image but not in the MSX 
E-Band could be due to the lower sensitivity of the MSX instru- 
ments. 

We include the main IR parameters in Table |4] The IR lumi- 
nosity hiR, the interstellar dust mass Md, and the dust temper- 
ature Td were derived assuming a distance D = 2.5 ± . 5 kpc . 
The IR luminosity was estimated following IChan & FichI i 1 9951) 
as Lir{Lq) = 1.58 5iR(Jy) (kpc), where Sir is the inte- 
grated flux given by Sir = 1.3(512 + S25) + 0.7(5'25 + Sm) + 
0.2(560 + 5ioo)- In order to derive the amount of radiatively 
heated dust, we have used the equation ( lDrainelll990h MdiM^) = 
m„56o ^^(kpc) (Bi'^^ - 1), where m„ = 0.30 x 10"® g 
cm~^ Hz~^ for the adopted value of n (n — 1.5) and _B„ is 
the modified Planck function, _B„ = 1.6 67^+" S im/Sm- As for 
the dust temperature, it was obtained as (lDraindfl99Ql Td{K) = 
(95.94/lnBn). 



3.3 The neutral hydrogen 

As mentioned in Section[T] molecular emission related to Sh2-I73 
was found at —34.5 kms^^ ( Blitz et al., 1982), which coincides 
with the Ha line velocity ob tained for the ionized gas, —34.3 ± 
0.3 kms"^ jFich et alj[l990l) . Thus, although the whole Hi cube 
was inspected looking for a structure related to the ionized region, 
special attention was paid to the velocity interval — 34±15 kms~^. 

Figure [SJj shows the mean brightness temperature image 
within the velocity range -30.3 to -25.4 km s-^ The presence of 
a small region of low Hi emission centered at (/, fc) = (1 19?5, — 1?0) 
is evident. This H I void is surrounded by a region of enhanced H I 
emission. All of the seven stars but one (ALS 6155) listed in Table 
[T]lie inside the H I minimum. ALS 6155 is seen projected onto the 
southern border of the cavity. In Figure [5}?, the 1420 MHz radio 
continuum image is overlaid in contours for comparison. The po- 
sitional coincidence is excellent, indicating that the H I feature and 
the Hll region Sh2-173 are quite likely physically related. 

A number of parameters characterizing this H I structure can 
be derived. Given its angular radius of about 14' , we calculate 
a physical size of 10 ± 2 pc. The expansion velocity was es- 
timated from the velocity range spanned by the Hi feature as 
(«2 — '^i)/2 + 1.3, where v\ and «2 are the minimum and max- 
imum velocities where the structure is observed, and the extra 1.3 
kms~^ takes into account that the Hi caps could be present at 
vi_ — 1.3 kms^^ and V2 + 1.3 kms^^. The systemic velocity of 
the structure corresponds to the velocity where the dimension of 
the cavity is largest. The mass of the Hi ring can be calculated 
under the assumption of optically thin H I emission. Following the 
procedure described by Pineault ( 1998), the Hi mass is given by 
Mhi(M0) = 1.3 X 10-3 Dgp^ At;J7„„2ATB , where Dg^^ 
is the distance to the H I structure and Aw the velocity width over 
which the H I cavity is detected, in km s"^. 

The meaning of ^am? ™d ATb depends on whether we want 
to find the missing mass in the cavity or the excess mass in the shell. 
In the first case ATb is the brightness temperature difference be- 
tween the void and the background, and ^am^ is the cavity solid 
angle in arcmin^ . In the second case ATb is the brightness temper- 
ature difference between the shell and the background, and r2a,„2 



Unveiling the birth and evolution ofSh2-l 73 5 

is the solid angle of the shell. All the derived parameters are listed 
in TableU 



3.4 CO emission distribution 

Figure |6] shows the "CO(l-O) emission distribution within the ve- 
locity interval from -27.4 to -40.6 kms^^ in steps of 3.3 kms"^ 
(contours) overimposed onto the DSS2-R optical image (grey- 
scale). Each panel is the average of four consecutive CO maps. 

In the velocity interval from -27.4 to -30.7 kms^^ , the CO 
emission partially borders the brightest optical emission arc at (/, b) 
~ (1 19? 25, -0?75). In the second velocity interval, from -30.7 to - 
34.0 km s~^, the CO structure extends towards the south increasing 
in intensity at lower galactic longitudes. Low molecular arc-shaped 
emission is distinguished surrounding the optical nebulosity at (/, 
b) ^ (I19?6, -0?75). The molecular emission distribution in the 
velocity range from -34.0 to -37.3 kms"^ changes in compar- 
ison with the previous images. Several different CO features are 
observed, being the dominant one a CO clump which spatially co- 
incides with the area lacking optical emission at (/, b) ~ ( 119^4, 
-0?85). The fact that radio continuum emission is detected in this 
region (see Fig. I2J3), reinforces the possibility that this molecular 
clump is the responsible for the observed optical absorption feature. 
In the last velocity interval, from -37.3 to -40.6 kms^^, the CO 
structures seen from -27.4 to -34.0 kms^^ disappear and most 
of the CO emission appears projected onto the optical absorption 
features situated in the central part of the H II region. In two of the 
velocity channels corresponding to this velocity interval, hardly de- 
tected (3 rms) molecular emission is seen projected onto (/, b) ~ 
(1I9?57, -1?0), in coincidence with the other area where optical 
emission is absent. This CO structure is indicated by a dotted con- 
tour in Fig.|6]. 

In Fig. |7] the averaged CO emission distribution in the veloc- 
ity interval from -27.4 to -39.0 kms~^ in comparison with the 
ionized gas, the dust, and the H I distributions in the environs of 
Sh2-173 is shown. The left panel of Fig. |7] shows the averaged 
CO emission distribution (contours) overimposed onto the DSS2-R 
optical image (grey-scale). There are two main clouds of molecu- 
lar gas probably related to Sh2-I73. Towards lower galactic lon- 
gitudes, molecular gas partially encircles the H II region (cloud A 
from here on). The dense molecular gas in this region has prob- 
ably slowed down the ionization front and the expansion of the 
Hll region in this direction. On the other hand, the area near (/, 
b) = (119°4, -0?85), which depicts relatively strong radio emission 
and is almost free of optical emission, coincides with a cloudlet of 
molecular material detected within the velocity range from -34.0 to 
^0.6 kms^^ (Cloud B in Fig.|7]l. This molecular gas and the as- 
sociated dust should be in front of the ionized material to originate 
the observed optical absorption. The velocity of cloud B (~ —38 
km s~^) is compatible with this material being pushed towards the 
observer. The difference in radial velocity between cloud A and 
cloud B allows us to get a rough estimation of the expansion ve- 
locity of the molecular gas, which results to be 4 ± 2 kms~^. The 
middle panel of Fig|7] shows an overlay of the image at 8.3 (xm 
(grey-scale) and the CO emission distribution (contours). The mid- 
infrared and molecular emissions present an excellent morpholog- 
ical correspondence. The molecular gas encircles the bright MSX 
rim. A fainter CO emission feature is observed towards the infrared 
emission detected at (/, b) ^ (119? 6, -0?75) at 8.3 fim. The right 
panel of Fig. [7] shows a comparison between the CO and H I emis- 
sion distributions, showing that the CO structure (cloud A) is lo- 
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Figure 3. HIRES IRAS maps of the same area as in Fig.|2] The small triangles indicate the position of the seven OB stars listed in Table[T] Grey-scale goes 
from 3 to 10 MJy sr^^ (dark grey) at 12 /xm, from 7 to 13 MJy sr~l at 25 /xm, from 17 to 55 MJy sr~l at 60 /im, and from 70 to 160 MJy sr"^ at 100 /xm. 



cated slightly further from the ionized gas as compared with the 
Hi shell. 

Summing up, the comparison between the spatial distribution 
of the ionized and neutral gas, together with the dust, reveals a re- 
markable morphological correspondence. The emission of the ion- 
ized gas is partially encircled by an arc shaped feature of mid- 
infrared emission and an outermost feature of molecular gas. This 
stratified distribution suggests that the ionization front is bordered 
by a PDR, detected as the bright rim at 8.3 ^m. 

Table|4]summarizes the physical parameters of the two molec- 
ular clouds shown in Fig. |7] (named cloud A and cloud B) as 



derived from the CO data. The masses were estimated by in- 
tegrating the CO line intensity as Wco = / T(CO)dv, where 
T(CO) is the average temperature of the molecular gas over 
the velocity interval where each cloud is observed. To calcu- 
late the H2 column density, the relation ship X = A^(H2)/Wco 
of 1.9 X 10^° cm~^(Kkms~^ )~^ jOrenier & Lebrunl 1 19901: 
iDigel. Hunter. & Mukherie^ll995h was considered. The molecu- 
lar mass was derived from M[Mq\ = 4.2 x 10'^° N{H2)D'^ A, 
where D is the distance in pc and A is the solid angle in steradians. 
Assuming a spherical geometry for both clouds, the volume density 
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Figure 4. Left panel: MSX Band A (8.3 ^m) image of the area around Sh2-173. Right panel: Overlay of the 8.3 fj,m image (grey-scale) and 1420 MHz 
emission (contours). Contour levels are at 6 and 7 K. The small triangles indicate the position of the seven stars listed in Table[T] 
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Figure 5. a) Averaged Hi distribution in the velocity interval —30.3 to —25.4 kms~^ showing the Hi minimum related to Sh2-173. Grey-scale goes from 
35 to 65 K. b) Overlay of the averaged H I distribution (grey-scale) and the 6 K radio continuum level at 1420 MHz showing the close coirespondence between 
the H II region and the H I cavity. The small triangles indicate the position of the stars listed in Table[T] 



can be estimated as n = 9.79 x M[Mq]/ , where R is the radius 
of the molecular cloud in pc. 

In order to obtain an estimate of the visual absorption Av 
produced by the gas in the direction of Sh2-173, we have used 
both the H I and CO data cubes. The column densities A'^(H l) and 
N{H2) were obtained by integrating the data cubes in the veloc- 
ity interval from to —40 kms^^ in a small box area (~ 2') 



over cloud A. The total hydrogen column density was obtained as 
Nh = A/^(Hl) + 2 N{H2)- Using the standard conversion factor, 
the total visual extinction was es timated as Av = 5.2 x 10~^^ Nh 
( iBohlin. Savage. & Drake 1978), which in this case yields Av = 
4 ± 1 mag. For comparison, the visual extinction towards the OB 
stars listed in Table[T] where n o molecular emissio n is observed, is 
in the range from 1.3 to 2 mag teusseil et al.ll2007l) . 
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Figure 6. Ionized and molecular gas in the environs of Sh2-173. The images display overlays of the CO emission distribution in the velocity range from -27.4 
to ^0.6 km {contours) and the DSS-R optical image (grey-scale). Contour levels are 0.5 K (= 5 cr), and from 1 to 5 K in steps of 1 K. Each panel shows 
the average of four consecutive CO maps. The dotted contour shown in the last velocity interval corresponds to the 0.6 K level of the image at ^0.2 km s~ ^ . 



4 STELLAR FORMATION ACTIVITY 

The multiwavelength study of Sh2-173 reveals that it is a classical 
expanding H II region containing gas and dust, partially encircled 
by a patchy molecular emission. The interface between the ion- 
ized and the molecular gas is clearly evident with the detection of 
a PDR. The observed morphology leads to the question of whether 
the expansion of the ionized region could have induced the forma- 
tion of a new stellar generation in the neighboring molecular enve- 
lope. In order to analyze this possibility, we have made a systematic 
search for young stellar object (YSO) candidates in the area. 

To look for primary tracers of stellar formation activity, we 
used the MSX6C Infrared Point Source Catalogue dSgan et al.l 
l2003h in Bands A (8.3 ^m), C (12.1 ^m), D (14.7 ^m), and 
E (21.3 nm); the 2MASS All-Sky Point Source Catalogue 



dCutri et al.ll2003h in bands J (1.25 nm), H (1.65 \im), and K 
(2.17 vim); and the IRAS Point Source Catalogue^. 

YSO candidates were searched in a 1° x 1° box area centered 
at (/, b) = (1 19?36, -0?9 ). A total of 30 IR AS point sources were 
found projected onto the analyzed region, ljunkes. Ftirst. & Reichl 
l ll992h 's conditions for YSOs are: Sioo ^ 20 Jy, 1.2 ^ ^ ^ 6.0, 

^ 1, and Qso+Qioo ^ 4, where Sa and Qa are the flux density 
and the quality of the IRAS flux in each of the observed bands, 
respectively. Only 2 of the 30 observed IRAS sources fulfill these 

conditions and may be classified as protostellar can didates. 

The MSX sources were classified based on the lLumsden et alj 

(l2002h's criteria. The selection of the sources is made taking into 
account their loci in the (F21/F8, F14/F12) diagram, where Fa 

1 1986 IRAS catalogue of Point Sources, Version 2.0 (11/125) 
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Figure 7. Dust and gas in the environs of Sli2-173. In the tliree panels the contours show the CO emission distribution within the velocity range from -27.4 to 
-39.0 kms^^. Contour levels ai'e from 0.2 to 1.8 in steps of 0.4 K. Left panel: The CO contours are superimposed on the DSS2-R optical image (grey-scale). 
Middle panel: The CO contours are superimposed on the MSX Band A image. Right panel: The CO contours are superimposed on the averaged H I distribution. 



denotes the flux in each band. According to these criteria, mas- 
sive young stellar object (MYSO) candidates have F2i/Fg > 2 
and Fi4/Ji2 > 1, while compact Hll regions (CHll) present 
F21/F8 > 2 and F14/F12 < 1. Evolved stars occupy the region 
where F2i/Fg ^ 2 and F14/F12 ^ 1. We have found 40 MSX 
sources projected onto the area, however neither of them satisfy 
the criteria to be classified as a MYSO nor a CH II. 

To look for tracers of stellar formation activity in the 
2MASS catalogue, we have adopted the criteria developed by 
iKerton et al 1 ( l2008h . Based on a sample of YSOs having low (< 
2 Mp)) and intermediate (2 < Mq < 9 ) masses extracted from 
iKenvon & Hartmannl ( Il995h (T Taur i) andlThe. de Winter. & PerezI 
( I1994 ') (Herbig/AeBe). respectivelv. lKerton et alj ^2008*) inferred a 
colour criteria to analyze the presence of stellar formation activ- 
ity in the surroundings of the H II region KR 140. They scaled the 
data of the sample to the distance of KR 140 (2.3 kpc) and the cor- 
responding visual extinction (Av = 5.5 mag). Using the 2MASS 
Point Source Catalogue, and taking into account different photo- 
metric qualities and the adopted colour criteria, these authors di- 
vided the YSO candidates in four groups, named PI, P1+, P2, and 
P3. PI and P1+ sources have the best photometric quality (ph.qual 
= A, B, C, or D). PI sources lie below the reddening vector asso- 
ciated with an 06 V star, while P1+ sources lie in the overlapping 
region of T Tauri and main sequence stars. The P2 group contains 
sources which have not been detected in the J-band (ph_qual = U 
in this band). In this case, the adopted J-magnitude is a lower limit 
and consequently the corresponding value of (J-H) is a lower limit 
as well. Thus, the actual position of the P2 sources in the (J-H)- 
axes is likely to be towards higher (J-H) values. Finally, the sources 
belonging to the P3 group have J and H magnitudes that are lower 
limits, so their colour (J-H) can not be considered. These sources 
can not be included in the CC diagram, and their real position on 
the CM diagram may be towards higher values of (H-K). To ana- 
lyze the presence of 2MASS YSO candidates in the surroundings 
of Sh2-173 we have adopted the Kerton et al.'s criteria scaled to the 
distance adopted for Sh2-173 (2.5 kpc) and the visual absorption 



Table 5. Selection criteria for 2MASS YSO candidates. 



Classification 




Colour criteria 


PI 


{J 


( J - H)> 0.835 
~ H)< 1.7 (//-if) -0.087 


P1+ 


(J 
{J- 


1.135 <{J -H)< 1.435 
~ H)> 1.7 (//-if) -0.087 
- //)< 1.7 (// - if) + 0.3685 
K > 14.5 


P2 


(J- 


- H)l < 1.7 (H - /f)- 0.087 
(H - K)> 0.91 


P3 




(H - K)l > 0.91 



The underscript "L" indicates a lower hmit. 



in this direction (Av = 4.0 mag). The modified criteria are listed in 
Table m 

A total of 23988, 135, and 134 sources from the 2MASS cat- 
alogue were examined in order to find PI or P1+, P2, and P3 
candidates, respectively. From the 23988 sources, only 81 and 10 
were classified as PI and PI +, respectively. Out of the 135 can- 
didates, 19 were identified as P2 sources, while only 5 out of the 
134 candidates were classified as P3 sources. Thus, a total of 115 
2MASS sources were found to be YSO candidates. However, 2 
of these sources had to be excluded from the sample since they 
were identified as part of the galactic contamination in the 2MASS 
catalogue (i.e. galcontam=2). Figure [8] shows the CC (top panel) 
and CM (bottom panel) diagrams for the YSO candidates pro- 
jected onto the molecular gas associated with Sh2-173. The CM 
diagram also shows the position of the Sh2-173's exciting stars. 
The redd ening vectors for an early type (09 V) and a late-type (MO 
111) stars ( lTokunag3l200(]h are represented by t wo parallel lines in 
the CC diagram using extinction values from IRieke & Lebofskvl 
( Il985l) . For the location of th e main sequence we h ave used for 
the 0-stars the values given bv lMartins & Pled ( |2006|) . and for the 
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Figure 9. Infrai'ed point sources classified as YSO candidates overlaid onto 
the CO distribution averaged between the velocity range from -27.4 to - 
49.8 kms~^. IRAS, PI, P1+, P2, and P3 sources are indicated as crosses, 
circles, squares, diamonds, and triangles, respectively. The locations of the 
four groups, gl, g2, g3, and g4, described in the text are indicated. 

late-type stars (B to M) t he ones given bv lTokunagal (l2000l) and 
iDrilling & Landolll j2000l) . The main sequence obtained for a dis- 
tance of 2.5 kpc is indicated by filled diamonds in the CM diagram. 
The location of the YSO candidates in the CC diagram i s in con- 
cordan ce with the numerical models developed by Robitai lTe et al.l 
( l2006h . 

Figure [9] shows the location of the YSO candidates with re- 
spect to the ^^CO emission distribution in the velocity range from 
-27.4 to -49.8 kms~\ IRAS, PI, P1+, P2, and P3 sources are 
indicated as crosses, squares, circles, diamonds, and plus signs, re- 
spectively. As can be inferred from the figure, 65 of the 115 YSO 
candidates appear projected onto molecular gas, from which 46 are 
located onto molecular clouds probably related to Sh2-173 (clouds 
A and B). There is a clear concentration of sources towards the 
brightest parts of the CO emission associated with Sh2-173. Based 
on their relative location, we have separated them into four groups: 
gl, g2, g3, and g4 (see Fig.[9}. The groups gl and g2 contain 8 and 
20 2MASS sources, respectively, and are located in the periphery 
of the optical nebula, onto the PDR. The group g3 is conformed by 
sources located towards lower galactic longitudes, onto cloud A. 
The sources of the group g4 are located onto cloud B. 

The main data of the IRAS and 2MASS sources that are seen 
projected onto molecular gas are listed in Table|6] It is worth men- 
tioning that an inspection of the CO data cube reveals that there is 
only one velocity component in the direction of Sh2-173. This fact 
strongly suggests that the YSO candidates suffer from extinction 
due to the molecular material related to Sh2-173. 

To sum up, based on different infrared colour criteria and on 
the location of the molecular emission, we have identified 46 YSO 
candidates probably associated with Sh2-173. 

The presence of a PDR is a clear proof of the interaction be- 
tween the H II region and the molecular gas, since it is a region 



directly affected by the expansion of the ionized gas and by the UV 
photons of the exciting stars. Thus, the relative location of the YSO 
candidates with respect to the PDR may be considered as an indi- 
cation of the role that the H II region is playing in the formation 
of the new stars. The probability for a source of being triggered by 
the action of the H II region is higher whe n it is located inside the 
region of direct influence of the H II region jKerton et al.l2008l , and 
references therein). 

Figure|9]shows that there are two groups of sources, named gl 
and g2, located adjacent to the PDR, while the group g3 is located 
further, towards lower galactic longitudes. The possibility for the 
group g4 of being located near or onto a PDR can not be discarted. 
A face-on PDR, which could not be detected, may be present be- 
tween the ionized region and the molegular gas. 



5 THE BIG PICTUIfE 

As mentioned in Section[T] Sh2-173 lies in the b order of a large H I 
shell observed in the Perseus arm Jpich 1986; Moor & Kissll2003h . 
Using CGPS data jTavlor et alj |2003). we have i nspect ed the area 
looking for the H I structure observed by IFichI ( Il986h using low 
resolution data. 

An inspection of the whole data cube reveals the presence of 
an Hi minimum which is first seen near —20 kms~^ and which 
can be followed up to —50 kms^^. Figure [Tol shows the Hi dis- 
tribution averaged in the velocity range from -20 to -50 kms^^, 
where the large Hi cavity is evident centred at (Z, b) = { 117? 8, 
+ 1?5), as well as the Hi shell encircling it. Following the stan- 
dard nomenclature for galactic shells, we designate this structure 
as GSH117.8-Fl.5-35. 

Along the mentioned velocity range, both the Hi minimum 
and the Hi shell are well defined as a large structure, although, 
as in many other H I shells, the typical behavior of an expanding 
structure (i.e. a ring that start as a cap, reach its greatest angular 
dimension at the systemic velocity and then shrink back to a cap) 
is not completely observed. In this case the receding cap is miss- 
ing, while the approaching cap is detected, at about -50 kms"'^. 
The absence of the receding cap could be due to confusion effects. 
Other possible explanation is that the receding part may be expand- 
ing into a lower density medium, as it would happen if for example 
it were located near the border of the Perseus arm. On the other 
hand, as expected when dealing with high resolution H I data, there 
are several small-scale structures inside and around the large ring. 
These smaller H I features change in shape and size along the veloc- 
ity range where the large shell is defined, and are possibly related 
to the individual action of every massive star present in the area. In 
this paper we will just focus on the large scale structure. 

The location of the stars belonging to Cas 0B5 
jGarmanv & Stencell 1 19921) are indicated in Fig. (TO] showing 
that most of them lie outside the large H I minimum. Particularly, 
two 0-stars belonging to CasOB5 are seen projected towards a 
smaller H I minimum centre d at (/, b) = (AlT.'i , -1?3). The position 
of the WN4 star (WR 159) ( lNegueruelj|2003h is indicated with a 
larger circle. 

In Figure[TT]a 1420 MHz radio continuum image of the same 
area displayed in Fig.[TO]is shown. The ellipse delineates the loca- 
tion of the H I shell. A minimum in the radio continuum emission in 
coincidence with the H I cavity is observed, although the presence 
of two SNR (CTB 1 and Gl 16.5-1-1.1) makes it less evident. Several 
Sharpless H II regions are projected onto the thi ck H I shell , being 
Sh2-173 one of them. Based on Ha observations, lFichetal.l l fT990l) 
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Table 6. YSO candidates projected onto the molecular clouds A and B. 



IRAS sources 



# 


Designation 


(lb) 




Huxes [Jy] 






Group 








12 jum 


25 //m 


60 //m 


100 //m 


For D = 2.5 kpc 




1 


00158+6128 


119? 057 -0? 87 


0.92 


2.29 


8.78 


23.29 


~ 180 


g3 


2 


00179+6121 


119? 282 -1?02 


0.32 


0.3 


5.89 


23.89 


~ 110 




2MASS sources 


# 


Designation 


(Lb) 


J 


H 


K 


(H-K) 


(J-H) 


Group 


PI 


3 


00 1 83077+6 145567 


119?047 -0?857 


14.317 


13.132 


12.041 


1.091 


1.185 


g3 


4 


00211376+6146520 


119?368-0?881 


16.892 


15.932 


15.317 


0.615 


0.96 


g4 


5 


00193852+6143342 


119? 175 -0?913 


16.816 


15.858 


15.07 


0.788 


0.958 


g2 


6 


00194881+6147427 


11 9? 204-0? 847 


15.799 


14.91 


14.322 


0.588 


0.889 


g2 


7 


00202395+6154406 


119?286-0?740 


17.227 


16.246 


15.499 


0.747 


0.981 


gl 


8 


00201448+6154455 


119?268-0?736 


16.76 


15.899 


15.007 


0.892 


0.861 


gl 


9 


00212578+6143577 


119?386-0?932 


16.827 


15.967 


15.356 


0.611 


0.86 


g4 


10 


00191976+6147310 


119? 146-0? 843 


16.041 


15.104 


14.472 


0.632 


0.937 


g2 


11 


00212969+6149057 


119?404-0?848 


17.11 


16.16 


15.359 


0.801 


0.95 


g4 


12 


00193126+6148555 


119? 172-0? 823 


16.706 


15.822 


14.851 


0.971 


0.884 


g2 


13 


00192242+6143509 


119? 144-0? 904 


16.382 


15.458 


14.68 


0.778 


0.924 


g2 


14 


00191965+6147518 


11 9? 147-0? 837 


16.709 


14.76 


13.56 


1.2 


1.949 


g2 


15 


00202211+6158243 


119?290-0?678 


16.455 


15.518 


14.651 


0.867 


0.937 


gl 


16 


00193625+6146178 


119? 176-0? 867 


16.28 


15.413 


14.767 


0.646 


0.867 


g2 


17 


00204518+6157266 


119?333-0?699 


14.873 


13.387 


12.373 


1.014 


1.486 


gl 


18 


00195258+6148460 


119?213-0?830 


16.51 


15.633 


14.967 


0.666 


0.877 


g2 


19 


00212565+6148318 


119?395-0?856 


16.559 


15.339 


14.464 


0.875 


1.22 


g4 


20 


00195062+6145397 


119?203-0?881 


15.993 


15.072 


14.361 


0.711 


0.921 


g2 


21 


00194259+6147191 


119?191 -0?852 


16.584 


15.694 


14.876 


0.818 


0.89 


g2 


22 


00192579+6146125 


11 9? 156-0? 866 


15.701 


14.751 


14.115 


0.636 


0.95 


g2 


23 


00213187+6149382 


119?409-0?839 


16.106 


14.707 


13.742 


0.965 


1.399 


g4 


24 


00212247+6148506 


119?389-0?850 


16.616 


15.265 


14.375 


0.89 


1.351 


g4 


25 


00204129+6159065 


119?329-0?671 


16.997 


16.016 


15.358 


0.658 


0.981 


gl 


26 


00193569+6148554 


119?180-0?824 


16.703 


15.597 


14.857 


0.74 


1.106 


g2 


27 


00191679+6142153 


119?130-0?929 


16.654 


15.789 


15.092 


0.697 


0.865 


g2 


28 


00202147+6155100 


119?282-0?731 


16.356 


15.48 


14.887 


0.593 


0.876 


gl 


P1+ 


29 


00210933+6146474 


119?359-0?881 


16.811 


15.655 


15.132 


0.523 


1.156 


g4 


30 


00192372+6147592 


119?155-0?836 


16.992 


15.674 


15.015 


0.659 


1.318 


g2 


31 


00211236+6149088 


119? 370-0? 843 


16.609 


15.287 


14.684 


0.603 


1.322 


g4 


P2 


32 


00182469+6148202 


119?041 -0?816 


>16.575 


16.038 


15.026 


1.012 


>0.537 


g3 


33 


00192419+6148555 


119? 158-0? 821 


>17.416 


15.415 


13.609 


1.805 


>2.001 


g2 


34 


00193839+6144163 


119?176-0?901 


> 16.996 


15.743 


14.791 


0.952 


> 1.253 


g2 


35 


00212409+6149210 


119? 393-0? 842 


>16.361 


16.017 


14.359 


1.657 


>0.344 


g4 


36 


00202552+6156054 


119?292-0?717 


> 16.886 


15.338 


14.119 


1.219 


> 1.549 


gl 


37 


00192337+6148577 


119? 156-0? 820 


> 17.844 


15.739 


14.21 


1.529 


>2.106 


g2 


38 


00212265+6149170 


119?390-0?843 


> 16.938 


15.862 


14.888 


0.973 


> 1.077 


g4 


39 


00204593+6156560 


11 9? 334-0? 708 


>16.118 


15.215 


13.961 


1.255 


>0.902 


gl 


40 


00212591+6149307 


11 9? 397-0? 840 


> 14.951 


15.553 


14.619 


0.934 
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Figure 8. CC and CM diagrams of the 2MASS YSO candidates projected onto clouds A and B. PI, P1+, P2, and P3 sources are indicated by circles, crosses, 
asterisks, and squares, respectively. Upper panel: {H — K, J — H) diagram. The positions of the derredened early-type main sequence and giant stars are 
shown. The reddening curves for MO III stars (upper line) and 09 V stars (lower line) are indicated. The colours (J — H) corresponding to P2 sourc es are 
lower limits. Lower panel: {H — K, K) diagram. The position of the main sequence )Tokunag3l2000l ; rDrimng & LandolfcoOOtlMartins & Plezll2003 ) with 
Av = 4 mag at a distance of 2.5 kpc is indicated by filled diamonds. The locations of the OB stars related to Sh2-173 are indicated by their corresponding 
names. The reddening curve for an 03 star is shown with a dashed line. The colours {H — K) con'esponding to P3 sources are lower limits. 



derived the radial velocity of galactic H II regions. From their work, 
we found out that most of the Hll regions shown in Fig.[TT]have a 
radial velocity between -32 and -55 kms^^, in coincidence with 
the radial velocity in which the large H I shell is detected. Only two 
regions, Sh2-171 and Sh2-167, have radial velocities outside this 
range, of abo ut —12 kms~'^ and —66 k ms~^ , respectively. As 
mentioned bv lCazzolato & Pineault|j2003h . Sh2-180, Sh2-181, and 
Sh2-18 2, may be associate d with CasOB7. Sh2-179 is a planetary 
nebula jKerber et al.ll2003h . Besides the H II regions, four SNR are 



seen projected in the area, Gl 14.3+0.3, Gl 16.5+1.1, CTB 1, and 
G120.2+01.4 (Tycho). 

The HIRES 60 /im image of the same area is shown in Fig. 
[T2I The large-scale stmcture is similar to the infrared loop GIRL 
Gl 17+00 shown by Kiss et al. ( 2004). The infrared emission shows 
a distribution similar to the radio continuum one. A minimum in 
the 60 fim emission is evident inside the large H I shell (indicated 
by the ellipse). One of the main differences between this image 
and the one shown in Fig.[TT]is the absence of SNRs at infrared 
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Figure 10. Averaged Hi emission distribution between t he vel ocity range from -50 to -20 kms^^. Tlie circles and crosses indicate tlie position of the stars 
belonging to Cas OB5 according to lGarmanv & Stencell jl992l) . The small circles coirespond to O-type stars, while crosses mark the location of B stars. The 
large circle indicates the position of WR 159. The ellipse delineates the location of the H I shell. 



wavelengths. As expected, the infrared emission associated with 
H II regions is clearly observed. 

Figure [T3]shows the averaged CO emission distribution in the 
velocity range from -50 to -20 km s ~ ^ . The CO data were obtained 
from the Five College Radio Astronomical Observatory (FCRAO) 
CO Survey. The image reveals several molecular clouds sparsely 
located at the shell's border. 

It is generally believed that expanding supershells may trig- 
ger star formation at their edges. As the shell expands, it be- 
comes gravitationally unstable, providing a physical mechanism 
to form clouds and induce star formation along its periphery (e.g. 
iMcCrav & Kafatosll 19871 : lElmegreeij 'l998) . There are some obser- 
vational evidence confir ming the importa nce of the shell's evolu- 
tion in creating new stars. IPatel et al J ( 119981 ) analyzed the origin and 
evolution of the Cepheus bubble and concluded that the Cep OB2 
association was probably triggered by the action of an earlier gener- 
ation of massive stars and suggested that a third generation of stars 
has recent l y form ed as a consequence of the evolution of the shell. 
lOev et alj 0553) analyzed the W3AV4 Galactic star forming com- 
plex and suggested that it consists of a hierarchical system of three 
generations: a supergiant shell triggered the formation of IC 1795 in 
W3, which in turn triggered the star formation in its surroundings. 
They confirmed this scen ario based on the age sequence showed for 
the different generations. I Amal & CortiH2007h found a new H I su- 
pershell (GS263-02-I-45) which has the OB-association Bochum 7 



located in its border. They concluded that the association may have 
been bom as a consequence of the evolution of GS263-02-I-45. 

In our case, having in mind the scenario previously described, 
i.e. a large expanding Hi shell containing several Hll regions in 
its edge, and where at least one of them, Sh2-173, shows evidence 
of star forming processes taking place in its surrounding molecular 
material (see Section |4]l, we wonder if this could be other case of 
a hierarchical system of three generations. If this were the case, 
it should be reflected by a corresponding age sequence. In what 
follows we will briefly discuss possible age estimates for both, the 
Hll region Sh2-173 and the large Hi shefl GSH 1 17. 8-1-1.5-35. 

Given that the main exciting star of Sh2-I73 is still visible, an 
upper limit for its age can be estimat ed. The main-sequence life - 
time for an 09 star is about 5.0 Myr jSchaerer & de Kotejri997h . 
On the other hand, a lower limit to the age can be estimated consid- 
ering the sound crossing time, R/vs, where Vs is the sound speed 
in the ionized gas (~ 15 kms^^) and R is the radius of the H II re- 
gion. For Sh2-I73 the corresponding lower limit is about 0.6 Myr. 
As a different approach, the age of the H II region can be infe rred 
using the simple model described bv lDvson & Williamsl ( Il98(l) for 
an expanding ionized region in a uniform medium. The expansion 
of an H II region is highly dependent on the density of the surround- 
ing gas. As a rough estimate, the original ambient density can be 
obtained by distributing the total mass related to the structure (ion- 
ized, neutral atomic, and molecular) over a sphere of 12 ± 2 pc 
in radius, which yields rio = 41 ± 25 cm~"^. Under these condi- 
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Figure 11. 1420 MHz radio continuum image from tlie CGPS survey. The ellipse delineates the H I shell shown in Fig llQI The location of the Sharpless H II 
regions and the SNRs present in the area are indicated by circles. 



tions, we can infer that the H II region has been expanding during 
0.1 - 1.0 Myr. This range for the dynamical age, together with the 
age limits estimated above, suggest that the age of Sh2-173 could 
be between 0.6 and 1.0 Myr. Using the same simple model the ex- 
pansion velocity of the Hll region can be inferred. We obtained 
T4xp = 6 ± 3 kms~^, which is in agreement with the expansion 
velocity estimated for the molecular gas (4 ± 2 km s^^). Thus, the 
ionized and molecular gas are expanding at the same velocity, as 
expected according to the models of Hosokawa & Inutsuka ( 2006) 
when the molecular material has been accumulated behind the ion- 
ization front. 

A rough estimate for the age of the large Hi shell, 
GSH 117.8+1.5-35 , can be obtained using a simple model to de- 
scribe the expansion of a shell created b y a continuous injection of 
mechanical energy IWeaveretalll 19771) . In this way, the dynam- 
ical age of the Hi structure can be estimated as tdyn = 0.6 J? 
(pc)/'K;xp(kms^^), where R is the radius of the shell and 14xp 
the expansion velocity. The expansion velocity can be obtained as 
Kxp = {v2 — tJi)/2 + 1.3 kms~^ (see Section [331 ). Given that 
the large Hi shell is observed over some 30 kms^^, an expan- 
sion velocity of about 16 kms"'^ is inferred. Thus, considering 
R — 135 ± 20 pc (for a distance of 2 kpc, see the discussion be- 
low) and Vcxp = 16 ± 2 kms^^, an age of 5 ± 1 Myr is obtained, 
which is large as compared with the age estimate obtained for Sh2- 
173. 

Concerning the origin of the large shell, iFici] ( Il986h sug- 
gested that the three supernova remnants located within the H I 



shell have actively enlarged it, which was presumably previously 
created by the joined action of several stellar winds and/or several 
supernova remnants. Based on the location of Cas 0B5 on the sky, 
|Mo6r & Kiss (2003) associated the expanding shell with the OB 
association. However, an energetic analysis has not been made in 
order to test this possibility. Are the massive stars (O-type stars) 
belonging to Cas 0B5 capable of creating such a shell?. The first 
condition that should be checked is if the Hi shell and CasOB5 
are at comparable distances from the Sun. The systemic velocity 
adopted for the la rge shell is —35 kms~^ , while according to 
lHuinDhrevsl (ll976l), CasOB5 has Vlsr = -39.5 kms"\ On the 
other hand, iGarmanv & Stencell i ll 99 A determined for Cas 0B5 a 
distance modulus of DM =11.5 mag, which yields a distance D 
= 2 kpc. As mentioned in Section [T] the presence of non-circular 
motions in the Perseus arm is important and should be taken into 
account. An inspe ction of the observed radial velocity field of the 
Galaxy derived bv lBrand & Blitzl ( 1 1993b gives for the adopted sys- 
temic velocity a distance of about 1.8-2 kpc. Thus, we will consider 
a distance of 2 kpc for Cas 0B5 and the H I shell. 

A second condition that should be checked is if the en- 
ergy budget of the massive stars is enough to create the large 
shell. The mechanical energy output of the stars should be com- 
pared with the kinetical energy of the shell. As can be seen 
in Fig. [H four O-type stars (BD-^62 2299 (08), BD-^61 2550 
(09.511), BD+61 2559 (09V), and HD 108 (05 f)) and one Wolf- 
Rayet star (WR159 (WN4)) are seen projected inside the large 
Hi minimum. All these massive stars are supposed to belong to 



Unveiling the birth and evolution ofSh2-l 73 15 




laa- 130" iia- ils" lu' 

Galactic Longitude 



Figure 12. 60 /im HIRES IRAS data of the same area as Fig[TT] The ellipse delineates the H I shell shown in FigfTol 



CasOB5 jGaimanv & Stencell [l992l : iNeguemelal l2003h . The to- 
tal mechanical energy released by the stars during a time t is 
-Eiu = (E0.5 M v'^)t, where M is the mass-loss rate and Vao 
the terminal velocit y of the winds. T he corresponding values were 
taken from .Howarth & Prinial jl989l) for the O-type stars (log M 
= -5.73, -6.76, -7.17, and -7.36; and = 2800, 2200, 2100, 
and 2000 kms"^ for the 5, 08, 09, and 09.5 stars, res pec- 
tively) and from j van d er Huch t. Cassinelli. & Williainj ( Il986h and 
IPrinia, Barlow. & HowartH (1199(f) for the WR star (log M = -4.5 
and Voo = 1900 kms"'^). We obtained that the total energy input 
due to the stellar winds is about = 1.3 x 10*^ i(yr) erg- 

The total kinetic energy stored in the shell is Ek = 
0.5 Afsh Kxp, where Meh is the total mass in the shell and Vcxp 
the expansion velocity. We estimated the expansion velocity of the 
shell to be 16 ± 2 kms~^. Following the procedures described in 
Sections [3.31 and [TA\ we obtained the Hi and molecular masses in 
the shell. Each mass was computed from the corresponding aver- 
aged image covering the velocity range from -20 to -50 kms^^. 
Adopting for GSH 117.8+1.5-35 a distance of 2 kpc, the total swept 
up mass is Mah = (8 ± 4) x 10^ Mp, w hich is consistent with the 
value estimated bv lMoor & Kisi ( |2003l) of 7.5 x W^Mq (for their 
adopted distance of 2.5 kpc, which implies a mass of 4.8 x 10^ M0 
for a distance of 2 kpc). Thus, the kinetic energy of the shell is 
Ek = (2 ± 1) x 10^^ erg. 

The theoretical model of IWeaver et alj ( 1 19771) estimates that 
the efficiency of conversion of mechanical wind energy E^ into 
kinetic energy of the shell Ek is up to 20 % for the energy con- 
serving model. This implies that the massive stars have to be act- 



ing during at least 8 Myr to create the shell. This value is larger 
than the lifetime of an O-type star and suggests that the five mas- 
sive stars can not alone generate such a shell and other source s 
of energy should be considered. As suggested by FichI ( Il986h . 
the supernova remnants located within the Hi shell may be ac- 
tively enlarging it. Out of the four SNRs located in the area, 
Gl 14.3+0.3 is projected outsi de the Hi shell and it is located 
at a d istance of about 700 pc ( lYar-Uvaniker. Uvaniker. & KothesI 
l2004h . Gl 14.3+0.3 and CTBl ha ve a distance estimate of about 
1.6 kpc (Yar-Uvaniker et al.ll2004l). while Tycho is located at about 
2.4 kpc ihse. Koo. & Tatematsul 20041) . Thus, assuming that there 
are three SNR possibly related to the shell and that each of them has 
injected 10^^ erg, then the massive stars have to be acting during at 
least 5 Myr, which is consistent with the dynamical age estimated 
for the large shell. It is important to note that given the large er- 
rors involved in these estimations, the possibility that other sources 
of energy (stellar winds and/or supernova explosions) would have 
contributed in the formation of the shell can not be discarted. For 
instance, it is highly probable that the progenitors of the three SNR 
have contributed in the creation of the shell. 

All the evidence shown above suggests that Sh2-173 is part 
of a hierarchical system of three generations. The members of 
Cas OB5 would have formed the large H I shell, whose expansion 
would have triggered the formation of the Hll region Sh2-173, 
which in turn is triggering new stars in the surrounding molecular 
material. The fact that Sh2-173 is one of several young H II regions 
lying at the edge of the H I shell (see Figll lb reinfor ces this con- 
clusion. For example, Sh2-170 has been analyzed bv iRoger et al.l 
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Figure 13. Averaged CO emission distribution in the velocity range from -50 to -20 km s ^ . The elhpse delineates the H I shell shown in Fig llQI 



( |2004) . who derived an age of (2.5 - 5) x 10^ yr, while Sh2-16 8 
is considered as a compact H 11 region bv lKothes & Kerto j Jiooj) . 
A detailed study of the H II regions shown in Fig. [TT] would be use- 
ful to better understand the importance that the expansion of the 
large H I shell had in the process of triggering the formation of new 
generations of massive stars. 



6 CONCLUSIONS 

The Hll region Sh2-173 has been analyzed in order to study the 
phenomena associated with the interaction of massive stars and the 
interstellar medium, and specifically, the process of triggering star 
formation. The analysis of the available radio continuum data have 
enabled us to estimate some parameters that characterize the H II 
region, like the ionized mass, the emission measure, and the elec- 
tron density. From the radio flux densities estimated at different 
wavelengths, the thermal nature of the source was confirmed by the 
estimation of the spectral index a = O.OiO.l, ^ v" . We noted 
that the massive stars present in the field can provide the necessary 
UV photon flux to keep the H II region ionized and heat the dust. 
As expected for Hll regions, Sh2-173 displays a strong correla- 
tion between the 1420 MHz brightness temperature and infrared 
emission. An inspection of the H I images has shown the existence 
of an H I cavity which presents a good morphological agreement 
with the ionized gas. The H I feature is observed within the veloc- 
ity range -24.6 and -31.6 kms^^ and has a systemic velocity of 
—27.0 ± 1.3 kms~^ , compatible with the velocity of the ionized 



gas as obtained from Hq observations. The distribution of the MSX 
emission at 8.3 /im clearly shows the presence of a PDR bordering 
the brightest part of the H II region. Molecular gas probably related 
to Sh2-173 is observed in the velocity range -27.4 to -4-0.6 km s~^. 
Particularly, the CO feature named cloud A is clearly observed en- 
circling both the H II region and the PDR. Applying different colour 
criteria, 46 YSO candidates have been identified in the vicinity of 
Sh2-173. Their positions in the CC and CM diagrams, as well as 
their correlation with molecular emission highly suggest the pres- 
ence of an active star formation process in the molecular envelope 
of Sh2-173. 



An analysis of the H I distribution in a larger area shows that 
Sh2-173 is located at the edge of a large shell, together with several 
other H II regions. This morphology suggests that Sh2-173 may be 
part of a hierarchical system of three generations. In this context, 
the age of the different structures were estimated. Based on evo- 
lutionary models of Hll regions, we estimated for Sh2-173 a dy- 
namical age of 0.6 - 1.0 Myr. As for the large Hi shell, we have 
obtained a dynamical age tdyn = 5 ± 1 Myr. Based on energetics 
considerations, we have concluded that the massive stars belonging 
to Cas 0B5 could not alone create the shell, being highly probable 
that at least three supernova explosions have contributed in its for- 
mation. The obtained age difference between Sh2-173 and the Hi 
shell, together with their relative location, lead us to the conclusion 
that Sh2-173 may have been created as a consequence of the action 
of a strong shock produced by the expansion of GSH 1 17. 8-1-1. 5-35 
into the surrounding molecular gas. The analysis of the ages and 
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distribution of the other H II regions observed in the area would be 
useful to confirm this scenario. 
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